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ABSTRACT 



Aims. The prompt optical emission contemporaneous with the y-rays from y-ray bursts (GRBs) carries important information on the central 
engine and explosion mechanism. We study the time lag between prompt optical emission and y-rays in GRB 990123 and GRB 041219a, 
which are the only two GRBs had been detected at optical wavelengths during the ascending burst phase. 

Methods. Assuming profiles of prompt optical light curves are the same as the prompt y-rays, we simulate optical light curves with different 
time lags and compare them with the observed optical flux. Then the best fit time lag and its error are determined by chi-squared values. 
Results. We find that time lags between prompt optical emission and y-rays in GRB host galaxy rest-frames are consistent in the two GRBs, 
which is 5 ~ 7 s for GRB 990123 and 1 ~ 5 s for GRB 041219a. This result is consistent with a common origin of prompt optical and y-ray 
emissions in the two GRBs. Based on synchrotron cooling model, we also derive the parameters for the two GRBs. 
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1. Introduction 

The prompt optical emission contemporaneous with the y-rays 
from y-ray bursts (GRBs) carries crucial information on the 
central engine and explosion mechanism. Internal shock mod- 
els and external shock models have been proposed to explain 
the prompt optical emission (see e.g. Meszaros & Rees 1999; 
Liang et al. 1999). However, the origin of the prompt optical 
emission remains an open issue. 

Prompt optical emission have been found for several GRBs 
during the brief durations of bursts (e.g. Akerlof et al. 1999; 
Vestrand et al. 2005; Rykoff et al. 2005; Klotz et al. 2006). For 
two of them, i.e. GRB 990123 (Akerlof et al. 1999; Kippen et 
al. 1999) and GRB 041219a (Vestrand et al. 2005; Barfhelmy, 
S., et al. 2004), early-time optical observations were carried 
out during their bursting phase, including both ascending and 
descending parts. Hence, it is possible to study the time lag be- 
tween prompt optical emission and y-rays in these two GRBs, 
which could provide important clues to the origin of prompt 
optical emission. 
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2. Data Analysis 

As shown in Figure 1, in each GRB, there are only three 
points with a positive optical detection (crosses) during the 
burst phase and the profile of optical light curve remains un- 
known. To avoid redundant parameters, we simply assume that 
the prompt optical light curve profile is exactly the same as 
the y-rays, but with a time lag to be determined. Then we 
simulate optical light curves with different time lags and com- 
pare them with the observed optical flux. The redshift of GRB 
990123 is 1.6 (Andersen et al. 1999; Kulkarni et al. 1999). The 
redshift of GRB 041219a is unknown. Barkov & Bisnovatyi- 
Kogan (2005) gave a redshift upper limit as z < 0.12 based on 
a model of a dust reradiation of IR afterglow in the envelope 
surrounding the GRB source. In this note, we assume z — 0.1 
for GRB 041219a. 

We consider the following four models for time lag be- 
tween optical and gamma-ray emissions: 

1 . No time lag between optical and gamma-ray emissions in 
GRB 041219a; 

2. No time lag between optical and gamma-ray emissions in 
both GRB 990123 and GRB 041219a; 

3. Time lags between optical and gamma-ray emissions in 
GRB 990123 and GRB 041219a are independent. 
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4. Time lags between optical and gamma-ray emissions in 
GRB 990123 and GRB 041219a are the same in their host 
galaxy rest-frames; 

Results of x 1 tests for the four models are given in Table 1 . 
For model 3, x 2 values of different time lags are shown in 
Figure 2, and simulated light curves with corresponding lags 
of minimum x 2 values are shown in Figure 1 . x 2 values of dif- 
ferent time lags for model 4 are shown in Figure 3. 

As shown in Table 1, model 2 is rejected. Model 4 and 3 are 
better than model 1, indicating that a common time lag between 
prompt optical emission and y-rays in both GRBs is quite pos- 
sible. In model 3 (Figure 2), for GRB 990123, the best-fitting 
time lag is 12.2 s, and there is another local minimum of x 2 val- 
ues around 18.2 s. The 12 s lag corresponds to moving the sec- 
ond pulse to the second optical point, while 18 s corresponds to 
moving the first pulse to the second optical point. Therefore the 
host galaxy rest-frame time lag should be 4.6 ~ 7.2 s. For GRB 
041219a, the best-fitting time lag is 3.0^' 3 s with Icr error, and 
the host galaxy rest-frame time lag will be 0.6 ~ 5.0 s. The^ 2 
value is very sensitive to the time lag in GRB 990123, thus in 
model 4 (Figure 3), the best-fitting lag value is dominated by 
the contribution of GRB 990123. 



3. Discussion and Conclusions 

It was reported in some previous studies that the behaviors and 
origins of the prompt optical emission in GRB 990123 and 
GRB 041219a are totally different: In GRB 990123, the op- 
tical emission was uncorrelated with the prompt gamma-rays, 
suggesting that the optical emission was generated by a re- 
verse external shock arising from the ejecta's collision with 
surrounding material; while in GRB 041219a, the optical emis- 
sion was correlated with the prompt gamma-rays, indicating 
optical emission was generated by internal shocks driven into 
the burst ejecta by variations of the inner engine (Vestrand et 
al. 2005). However, our results show that time correlation be- 
tween prompt optical and gamma-ray emission are quite con- 
sistent in both GRBs. If GRB 041219a is a low-redshift object, 
rest-frame time lags of prompt optical emission behind y-rays 
in the two GRBs is similar, indicating a common mechanism 
to produce the prompt optical emission in the two GRBs. 

If both prompt optical and y-ray emissions come from 
cooling e ± pairs in the moving ejecta with a Lorentz factor 
T through synchrotron radiation, assuming velocities of e ± 
are isotropic, the frequency of emitted photon is around the 
peak frequency of synchrotron photons, which in the observer's 
frame is 



v « v m * 2 x 10 6 y 2 fiI7(l + Z ) Hz 



(1) 



where y is the Lorentz factor of e ± and B is the magnetic field 
strength in Gauss. 

The power radiated per one electron or positron in the ob- 
servation frame is 



The time lag between prompt optical emission and y-rays 
should be the typical lifetime of e ± which emit optical photons 



Has ~ „ 



3.3 x 10 8 (1 +z) 
yB 2 T 



s, 



(3) 



Assuming F = 300, which is a typical value in GRBs (see 
e.g. Liang et al. 1999), we can derive B, y and the total e ± 
number N e based on synchrotron cooling, as listed in Table 2. 
Further assuming an emission radius R = 2 x 10 16 cm (see e.g. 
Liang et al. 1999; Li et al. 2003), we can derive the synchrotron 
self-absorption frequency v a (Li et al. 2004) and the e ± column 
density S in the observer's frame 



v„ * 1 x 10 I6 L 2 { 7 r^7 4/7 fi! /7 /(l +z)Hz, 
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and 



I = 



N„ 



(4) 



(5) 



P = 1.1 X 10- 15 yVr 2 /(l + z) 1 erg/s. 



(2) 



As shown in Table 2, the self-absorption frequencies are 
less than the observed frequencies, thus prompt emissions can 
be observed. When absorption is negligible, the optical to y- 
ray flux ratio from synchrotron cooling with constant F and B 
should be F 5m k /F m keV = v 5000 A /v 100 kev = 2.5 x 10~ 5 . This 
is consistent with the observed value of GRB 041219a, which 
is 1.2 x 10~ 5 when there is zero lag between prompt optical 
emission and y-rays (Vestrand et al. 2005) and about 1.4x 10~ 5 
for our best-fitting model. However, in GRB 990123, the peak 
optical/peak y-ray flux ratio is 3.3 x 10~ 4 (Kulkarni et al. 1999; 
Akerlof et al. 1999) and the flux ratio changes to be 7.1 x 10 4 
for our best-fitting model when the time lag is 12.2 s, which is 
one magnitude larger than the synchrotron-cooling model pre- 
dicted value, indicating that such a simplified non-absorption 
synchrotron cooling model alone could not account for all the 
observed properties of GRB 990123. It is possible that an- 
other mechanism may be operating simultaneously which can 
modify the flux ratio. For example, down-Comptonization of 
gamma-ray photons in intervening electron clouds may be able 
to reduce the flux ratio substantially if some gamma-ray pho- 
tons are converted into optical photons which are also delayed 
by several seconds with respect to the gamma-ray emission, as 
proposed independently by Zheng et al. (2006). In this model, 
different covering factors of high density region around the cen- 
tral engine could lead to different optical to gamma-ray flux 
ratios. 
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Table 1. Results ofx 2 tests for different optical lag models 



model " 


x 1 


V 


POx 1 )" 


Lag 99 c (s) 


Lag 04 " (s) 


1 


7.22 


4 


0.12 







2 


5150 


6 


< io- 20 








3 


6.36 


4 


0.17 


12.2 


3.0 


4 


7.17 


5 


0.21 


12.2 


5.2 



" Refer to text for details. 

* The probability that a single observation from a x 2 distribution 
with degrees of freedom v will be larger than the x 1 values in column 
(2). 

c Time lag for GRB 990123. 
d Time lag for GRB 041219a. 

Table 2. Parameters of GRB 990123 and GRB 041219a based 
on synchrotron cooling model 



Parameter 


GRB 990123 


GRB 041219a 


Ref. " 


Obervation 








Redshift 


1.6 


0.1 


1,2,3 


Time lag (s) 


12.2 


3.0 




Lpeak.y (erg/s) 


1.5 x 10 53 


1.2 x 10 52 


2,4 


Lpeakfipt (erg/s) 


5 x 10 49 


1.5 x 10 47 


5,6 


Vl00 keV (Hz) 


2.4 x 10 19 


2.4 x 10" 




V 5000A (Hz) 


6x 10 14 


6x 10 14 




Assumption 








Bulk Lorentz factor 


300 


300 


7 


Radius (cm) 


2x 10 16 


2x 10 16 


7, 8 


Derived parameters 








B (Gs) 


28 


53 




Jy 


6.1 x 10 4 


2.9 x 10 4 




Jopt 


310 


140 




N e b 


3.6 x 10 51 


6.3 x 10 49 




V a ,y (HZ) 


6.5 x 10 14 


8.2 x 10 14 


9 


v appt (Hz) 


6.6 x 10 13 


3.2 x 10 13 


9 


X (cm" 2 ) 


2.6 x 10 23 


4.5 x 10 21 





" References: 1. Andersen et al. 1999; 2. Kulkarni et al. 1999; 
3. Barkov & Bisnovatyi-Kogan 2005; 4. Barthelmy et al. 2004; 5. 
Akerlof et al. 1999; 6. Vestrand et al. 2004; 7. Liang et al. 1999; 8. 
Li et al. 2003; 9. Li et al. 2004. 
b Calculated from L„ 
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Fig. 1. Observed and simulated optical and y-ray light curves. 
The upper panel is for GRB 990123 with time indicated in 
the top axis, and the lower panel is for GRB 041219a with 
time indicated in the bottom axis. The optical flux measured by 
Akerlof et al. (GRB 990123) or Vestrand et al. (GRB 041219a) 
are indicated by crosses, where error bars for detections are 
given as l<x values and non-detections are plotted as 2<x up- 
per limits. The corresponding best-fitting simulated values for 
model 3 described in the paper are plotted as circles. The y-ray 
light curves measured by the BASTE (GRB 990123) or Swift 
BAT (GRB 041219a) are plotted as solid lines. The best-fitting 
simulated optical light curves for model 3 are plotted as dashed 
lines. 
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Fig. 2.x 1 value of simulated optical flux as a function time lag 
for model 3 described in the paper. The upper panel is for GRB 
990123 and the lower panel is for GRB 041219a. The solid, 
dashed and dash-dotted horizontal lines indicate l<x, 2<x and 
3<x errors, respectively. 
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Fig. 3. x 1 value of simulated optical flux as a function time lag 
for both GRB 990123 and GRB 041219a in their host galaxy 
rest-frames, for model 4 described in the paper. The solid, 
dashed and dash-dotted horizontal lines indicate lcr, 2cr and 
3cr errors, respectively. 



